Interferometric Observation of the Highly Polarized SiO Maser 
Emission from the v = 1,J = 5 — 4 Transition Associated with VY 

Canis Majoris 

Hiroko Shinnaga^'^, James M. Moran^, Ken H. Young^, and Paul T.P. Ho^ 

hshinnagaScf a . harvard . edu 

ABSTRACT 

We used the Submillimeter Array to image the SiO maser emission in the 
f = l, J = 5 — 4 transition associated with the pecuhar red supergiant VY Canis 
Majoris. We identified seven maser components and measured their relative 
positions and linear polarization properties. Five of the maser components are 
coincident to within about 150 mas (~ 200 AU at the distance of 1.5 kpc); most 
of them may originate in the circumstellar envelope at a radius of about 50 mas 
from the star along with the SiO masers in the lowest rotational transitions. 
Our measurements show that two of the maser components may be offset from 
the inner stellar envelope (at the Sa level of significance) and may be part of 
a larger bipolar outfiow associated with VY CMa identified by Shinnaga et al. 
The strongest maser feature at a velocity of 35.9 kms~^ has a 60 percent linear 
polarization, and its polarization direction is aligned with the bipolar axis. Such 
a high degree of polarization suggests that maser inversion is due to radiative 
pumping. Five of the other maser features have significant linear polarization. 

Subject headings: masers — polarization — techniques: interferometric — stars: 
late-type — stars: individual (VY Canis Majoris) — radio lines: general 

1. INTRODUCTION 

The unusual star VY Canis Majoris (VY CMa), is located at the edge of the very large 
(~ 5° in diameter) Sharpless 310 HII region at a distance of 1.5 kpc (e.g. Herbig (1970); 
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Lada & Reid (1978)), and is one of the most intrinsically luminous stars known in the 
Galaxy with ~ (2 — 5) x IO^Lq (e.g. Le Sidaner & Le Bertre (1996)). Despite its high 
luminosity, it has a low effective temperature of ~ 2800 K, making it a red supcrgiant 
of spectral class M5Ib with a mass of ~ 25 Mq. Its unusually high mass loss rate for a red 
supergiant, M^, ~ (1 — 3) x 1O~^M0 yr~^ (Jura & Kleinmann 1990), is consistent with its 
thick circumstellar envelope and suggests that VY CMa may be close to the end of its life 
cycle. High resolution images taken at optical and infrared wavelengths reveal complicated 
circumstellar dust structures (Monnier et al. 1999; Smith et al. 2001). However, it is difficult 
to investigate the innermost circumstellar structures via the continuum emission because of 
the high opacity. 

VY CMa's extended outflowing envelope is traced by various species of molecular masers 
that conform to the standard paradigm for late type and supergiant stars whereby various 
species trace specific layers of the envelope for reasons involving density and excitation 
process (e.g., Moran et al. (1977); Rcid & Muhleman (1978); Marvel, Diamond, & Kemball 
(1998); Shinnaga et al. (2003); Miyoshi (2003)). One of the unusual characteristics of the 
emission from the circumstellar envelope is the high degree of linear polarization in the SiO 
V = lines. The only other stars that are known to show polarization in the SiO ground 
state are R Gas (Shinnaga, Tsuboi, & Kasuga 2002) and Orion KL/IRG2 source I (Tsuboi 
et al. 1996). 

We report the first imaging for the v = 1, J = 5 — 4 transition with high angular 
resolution and also describe the spectral characteristics of the strong linear polarization 
detected in this line. 



2. OBSERVATIONS AND DATA REDUCTION 

Observations of the SiO emissions in the v = 1, J = 5 — 4 transition (215.595950 GHz) 
towards VY GMa were made with the partially completed SubmilUmeter Array^ (SMA; 
Ho, Moran, & Lo (2004)), located on Mauna Kca, on 7 December 2002 under fair weather 
conditions (the zenith opacity r225GHz ~ 0.15 — 0.2). The phase tracking center was set 
at RA = 07'^ 22"* 58.^27, Deck = -25° 46' 03."4 (J2000) and the field of view (FOV) was 
57". A compact configuration (baselines ranging from 6.3 to 18.1 kA) with three antennas 
was employed for imaging the maser fine. The beam size was 13.9" x 5.6". The spectral 
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window was centered at 10.0 kms""^ in Vlsr- The frequency resolution was ~ 400 kHz. 
The typical double sideband system temperature was 500 K. We observed VY CMa along 
with a relatively bright nearby quasar, 0727-115 (~ 14° away), for the relative amplitude 
and astrometeic phase calibration, and Callisto for the absolute flux calibration. Precise 
phase calibration was achieved through the use of the maser component at 35.9 kms~^ as 
a phase reference. The passband calibration was done on Jupiter. Visibility data were 
calculated with the MIR IDL package (Scoville et al. 1993). The quasar amphtude showed 
no systematic variation with random fluctuation less than ~ 10% of the mean value, which 
reflects the accuracy of the amplitude calibration. 

Linear polarization measurements were accomplished by taking advantage of the diurnal 
rotation of the sky with respect to the axis of polarization sensitivity of the Array. The 
SMA antennas have alt-azimuth mounts, and the receivers, mounted at the Nasmyth focus 
following reflections from five intermediate mirrors, have fixed linearly polarized feeds. The 
final reflection before the receiver is achieved with a wire grid, which ensures that the overall 
polarization purity is better than 20 db. Because of the Nasmyth optics, the plane of the 
receiver polarization outside the antenna with respect to the local vertical is 45 — EL, where 
EL is the elevation angle of the source. The feed position angle on the celestial frame is 
(e.g.. Smart (1962)) 

PA(deg) ^45-EL + 3i /"K'^LAT) • MHA) 

cos[EL) 

where (f) is the latitude of the Array (~ 19°49'27") and HA is the source hour angle. In our 
experiment, HA ranged from —3 hours to 2.5 hours and PA ranged from —25 to 60°, which 
was sufficient to determine the linear polarization characteristics of the maser emission. 



3. ANALYSIS, RESULTS, AND DISCUSSION 

3.1. Line profile and linear polarization of the maser 

Figure 1 shows an integrated spectrum (vector average of all spectra) of the SiO high J 
maser (■?; = 1, J = 5 — 4). The noise level in the spectrum was 1.1 Jy. Seven major velocity 
components have been identifled by the number in the flgure. In addition to these, we also 
detected two weak high velocity components at the red shifted edge, at Vlsr ~ 48 and ~ 
56 kms~^. The brightest maser, component 2, has a very narrow line width (< 0.5 kms~^ 
at full width at half maximum) corresponding to a thermal Doppler width of about 600K if 
the maser is saturated. The line profile of the high J transition is quite different from those 
of the lower J transitions (e.g., Shinnaga, Tsuboi, & Kasuga (1999)), indicating different 
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regions for their origins. 

Figure 2 the flux densities of the seven major components at each PA. To determine 
the polarization of the masers, we fit the fiux density of each maser component, Sk, as a 
function of PA to the equation 

Sk^So,+Sp^cos[2{PA-i;k)] (2) 

where Sq^ is the average fiux density of maser component k, Sp^, is the polarized fiux density, 
and ipk is the position angle of the linear polarization. Table 1 shows the results of this 
analysis. The fractional linear polarization, listed in the rightmost column of Table 1, is 
given by iSp^/iSoj.. The fractional polarization percentages ranged from less than 10 percent 
(component 5) to about 60 percent (components 1 and 2). These values of fractional polar- 
ization are lower limits to the total polarization, since our measurement technique was not 
sensitive to circular polarization. 

Since component 5 showed no measurable change in flux density over the 5.5 hour dura- 
tion of the observations, we are confldent that the gain of the SMA was correctly calibrated 
throughout the experiment by the normal calibration on the nearby quasar, and that this 
component has no measurable linear polarization. To check this conclusion, and possibly 
improve the gain calibration, we performed the following analysis. After the linear analysis 
(based on Eq. (2)) to determine the three polarization parameters for each maser compo- 
nent, we examined the residuals at each time (PA). We estimated a gain factor correction 
for each time (six parameters) to minimize the overall deviations. We then solved again for 
the polarization parameters, and iterated the process several times. The solution converged 
quickly and stably, yielding the six gain factors: 1.00, 1.00, 1.04, 0.94, 0.99, and 1.03. The 
flnal polarization parameters were not signiflcantly different from those obtained in the flrst 
iteration (i.e., gain factors set to unity). The iterative solution reduced the overall rms noise 
level by about 20 percent and reduced the to essentially unity. This demonstrates that 
the measurement of the polarization parameters are accurate to within the quoted errors 
and are not subject to systematic errors due to calibration errors in the gain of the Array. 
It also shows that the a priori gain calibration based on the interleaved quasar observations 
was quite good, and validates the implicit assumption that the maser components did not 
vary intrinsically by more than a few percent over the period of observations. 

We tested the hypothesis that the polarization angles of the maser spots are the same. 
Under this assumption, the weighted mean polarization angle is 72 ± 4°. The reduced 
is 1.13, which strongly suggests that the polarization angles are the same and that the 
deviations from the mean of the position angles are not signiflcant, except possibly for 
components 3 and 6, which each deviate by about 1.5 a. 
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3.2. Angular distribution of the maser spots 

We assume that the individual masers components are unresolved by the SMA because 
their intrinsic sizes are less than ten milliarcsecond (e.g., see Reid & Moran (1986) for 
discussion of maser angular sizes). We tried to measure the angular distribution of maser 
spots with respect to the strongest maser component 2. To estimate the positions, we made 
an image for each spectral channel and fit it with a two dimensional Gaussian function, 
using AIPS tasks IMFIT and JMFIT. To improve the component position error estimates, 
we averaged the results over several channels that contribTitcd to the maser emission of 
each spot. Note that, for a point source, the determination of the peak of the image is 
statistically equivalent to a fitting of the phases (e.g. Thompson, Moran, & Swenson 
(2001)). In particular, since the analysis is based on the assumption that the individual 
maser components are point sources, the fact that the synthesized beam had high sidelobes 
was not important. Note that the visibility amplitude variations due to polarization do not 
bias the position estimate. 

Since we referenced the phase of each maser to the phase of the reference feature, the 
effects of any systematic errors were eliminated to first order. Any random phase errors 
between channels were reflected in the formal position errors. However, any systematic 
phase drifts between channels, caused by delay changes or electronic changes in the filter 
responses could have affected the measured positions. We have analyzed the phase versus 
time data for the reference quasar, spectrally resolved into four channels, and we found no 
evidence for any relative phase changes at the limit of about two degrees, which was limited 
by the measurement noise. The electronic stability and delay calibration of the system 
suggest that the relative phase changes are much smaller than this amount. A two degree 
systematic error would produce an error of less than 0.02" in relative position. 

Table 1 lists the positions of the masers, the relevant velocity ranges, and the mea- 
surement errors. The errors are close to the expected levels predicted by the formula 
59 = 0.5BW/snr (e.g. Kogan (1996)), where BW is the beamwidth in direction of RA 
or Deck, and snr is the signal-to- noise ratio. The snr of each component is its flux density 
divided by the measured noise of 1.10 Jy obtained from the integrated spectrum. 

Masers 3, 4, 5, and 7 are coincident with maser 2 within the 1-a error bars (150 mas in 
each coordinate or 200 AU). These masers are most likely associated with the SiO masers 
in the circumstellar shell at radius of 50 AU. Higher resolution measurements are needed 
to determine their exact positions in the shell. Masers 1 and 6 are offset in declination by 
1300 mas {2.5a) and 500 mas (4(t), respectively. These offsets might be real, but clearly 
have only marginal significance. Although there appear to be no J = 1 — transition masers 
beyond the 50 mas shell (Miyoshi 2003), it is conceivable that masers could exist in the 
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high J transitions in the bipolar outflow due to shock excitation. Masers 1 and 6 could be 
manifestations of such shocks. 

The positions of the maser spots are plotted in Figure 3, with the polarization indicated 
on each spot. We have placed component 2 between the two outflow lobes at the position set 
by the central star. If the assumption is correct, the uncertainty of the relative alignment of 
the two images is < 0.1". The relative alignment of the v — image and our maser image 
is uncertain to < 0.5", based on the positions of the reference quasars observed in these two 
measurements. Note that the position angle of the polarization of maser component 2 (71 
± 5°) and that of the mean of all the components (72 ± 4°) are nearly aligned with the 
outflow. This correlation suggests a relationship between the circumstellar envelope and the 
more extended bipolar outflow. Note that since the components are on both sides of the 
systemic velocity, they arise in both the approaching and receding sides of the envelopes. 

3.3. Continuum measurement 

We imaged the continuum emission in the line free channels. Because of the poor 
beamshapc of the Array, we could not determine the spatial structure of the continuum. 
However, we measured the total flux density to be 270 ± 40 mJy. The flux at 301 GHz and 
658 GHz measured with the SMA recently were 340 ± 10 mJy and 7.8 ±2.6 Jy, respectively. 
The SED obtained by combining our results with other measurements in the submillimeter 
band (Marshall, Leahy, & Kwok 1992; Sopka et al. 1985) give a spectral index for flux density 
of about — 3.0 for the frequency range 200 — 750 GHz. The emission is clearly almost entirely 
due to dust, with a small contribution from the star. The spectral index is slightly less ( — 
2.5) over the range 1,000 - 15,000 GHz (300 - 20 microns) (Le Sidaner & Le Bertre 1996). 

4. Pumping Mechanism of the SiO High J Maser 

Two theories have been proposed to explain the SiO maser pumping mechanism, es- 
pecially for the low J transitions: radiative pumping (Kwan & Scoville 1974; Deguchi & 
Iguchi 1976; Bujarrabal 1994a) and collisional pumping (Elitzur 1980; Doel et al. 1995). The 
positional coincidence of emission among different v > 1 vibrational states as observed with 
KNIFE (Miyoshi et al. 1994) seems to support the collisional pumping theory. 

On the other hand, the 60% linear polarization that we observe suggests radiative pump- 
ing (Nedoluha & Watson 1990). A fractional polarization this high for the J = 5—4 transition 
apparently requires anisotropy in the angular distribution of the IR radiation involved in de- 
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termining the populations of the masing states (Western & Watson 1984) .This anisotropy 
could occur, in principle, for either coUisional (through the escaping IR radiation) or radia- 
tive pumping. However, it occurs more naturally for radiative pumping since the incident IR 
radiation comes from the direction of the star, and is a larger effect in this case (Western & 
Watson 1983). Note that most theoretical analyses predict lower polarization in the higher 
J transitions, contrary to observations. 

The polarization direction of the radiation would be parallel or perpendicular to the 
magnetic field if gO, > R, F, where gQ is the Zeeman rate, 1.5 rad s^^ mG^^, F is the maser 
level decay rate, about 5 s~^, and R is stimulated emission rate, 4 x 10^® Tb Qm, where 
Th is the brightness temperature and Qm is the maser beam angle (Goldreich, Keeley, & 
Kwan 1973). Otherwise, the polarization would be aligned with the direction of the pump 
anisotropy. The VLB A image of the SiO masers in the v — 1,J— 1 — transition suggest 
that the spot sizes are about 2 mas (about 2 AU) (Miyoshi 2003), so the maximum brightness 
temperature is about 4 x lO^K. Taking Qm = 10~^ (e.g., Reid & Moran (1986)), the masers 
satisfy the gQ R criterion if i? > 10 niG. Since the magnetic field in the water masers in 
VY CMa has been estimated to be about 150 mG, based on a Zeeman interpretation of the 
observed circular polarization (Vlemmings, Diamond, & van Langevelde 2002), polarization 
would be expected to trace the magnetic field. Note that well ordered magnetic fields have 
been seen in other SiO masers (e.g. VY CMa: Shinnaga et al. (2003), TX Cam: Kemball & 
Diamond (1997)). 

We thank all SMA team members for their enthusiastic work and the Hawaiian people 
for allowing us to use their sacred mountain, Mauna Kea. We are grateful to Raymond 
Blundell, Shuji Deguchi, Ray Furuya, Liz Humphreys, Makoto Miyoshi, Scott Paine, and 
WiUiam Watson for helpful discussions. 
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Table 1. Parameters of Seven Maser Components 



Number 


Velocity'' 


Velocity'' 


Peak Intensity 


X Offsef^ 


Y OSsef" 


PA 


Degree of Linear 




(kms-i) 


Range 


(Jy) 


(arcsec) 


(arcsec) 


(degrees) 


Polarization (%) 


1 


42.1 


41.6 - 42.7 


7.4 


-0.23 ± 0.30 


-1.28 ± 0.50 


92 ± 30 


63 ± 31 


2 


35.9 


33.1 - 37.6 


63.9 


0.00 


0.00 


71 ± 5 


56 ± 4 


3 


26.9 


25.2 - 30.3 


33.6 


-0.01 ± 0.06 


-0.14 ± 0.09 


31 ± 24 


23 ± 6 


4 


24.1 


22.4 - 24.1 


23.7 


-0.02 ± 0.09 


-0.08 ± 0.13 


62 ± 21 


16 ± 4 


5 


16.1 


13.9 - 17.3 


23.0 


0.01 ± 0.10 


0.05 ± 0.15 




<10 


6 


12.2 


11.1 - 13.3 


20.2 


0.14 ± 0.07 


0.46 ± 0.10 


89 ± 12 


32 ± 6 


7 


0.9 


-0.3 - 2.0 


15.9 


-O.Olit 0.15 


0.13 ± 0.21 


78 ± 60 


29 ± 25 



*SR velocity of maser component peak (used for polarization analysis). 
''Range of velocities used for position analysis. 

•^The relative positions with respect to the brightest maser component, number 2. 
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Fig. 1. — An integrated spectrum of the SiO high J maser (v — 1, J = 5 — 4) associated 
with VY CMa made from the sum of all observations at various position angles. The velocity 
axis (radio definition) is measured with respect to the Local Standard of Rest. The stellar 
velocity is about 20 kms^^. There are probably additional components at 48 and 56 kms^^, 
which were too weak for detailed analysis. The numbers with arrows identify the seven major 
components. 
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Fig. 2. — The apparent flux density of seven velocity components as a function of the position 
angle of the feed polarization. The flux density of component 5 was almost constant over 
the PAs, the flux of the other six components changes as a function of PA because of their 
linear polarization. No fine gain adjustments were made in the data presented here (see 
Section 3.1). 
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Fig. 3. — Distribution of the SiO {v = 1, J = 5 — 4) maser components, overlaid on the VLA 
image of the SiO (f = 0, J = 1 — 0) emission (Shinnaga et aL 2003). The absolute position 
of maser 2 is RA(J2000) = T'^ 22™ 58.31" ± 0.03^ Decl.(J2000) = -25°46'3."2 ± 0."5. The 
area of each spot is proportional to its flux density. The line lengths are proportional to the 
polarized flux density. Only maser 6 (the most northerly component), and possibly maser 
1 (the most southerly component) have statistically significant offsets from the position of 
maser 2 (see Table 1). The mean polarization angle is 72 ± 4°. None of the deviations from 
this mean are significant, except perhaps for masers 3 and 6 (see Table 1). 



